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We present local extrema studies of two models that introduce a preferred direction 
into the observed cosmic microwave background (CMB) temperature field. In par- 
ticular, we make a frequentist comparison of the one- and two-point statistics for 
the dipole modulation and ACW models with data from the five-year Wilkinson Mi- 
crowave Anisotropy Probe ( WMAP). This analysis is motivated by previously revealed 
anomalies in the WMAP data, and particularly the difference in the statistical nature 
of the temperature anisotropics when analysed in hemispherical partitions. 

The analysis of the one-point statistics indicates that the previously determined 
hemispherical variance difficulties can be apparently overcome by a dipole modulation 
field, but new inconsistencies arise if the mean and the ^-dependence of the statistics 
are considered. The two-point correlation functions of the local extrema, £tt (the 
temperature pair product) and £pp (point-point spatial pair-count), demonstrate that 
the impact of such a modulation is to over-'asymmetrise' the temperature field on 
smaller scales than the wave-length of the dipole or quadrupole, and this is disfavored 
by the observed data. The results from the ACW model predictions, however, are 
consistent with the standard isotropic hypothesis. The two-point analysis confirms 
■ that the impact of this type of violation of isotropy on the temperature extrema 

statistics is relatively weak. 

From this work, we conclude that a model with more spatial structure than the 
dipole modulated or rotational-invariance breaking models are required to fully explain 
the observed large-scale anomalies in the WMAP data. 

Key words: methods: data analysis - cosmic microwave background. 
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1 INTRODUCTION 



Local extrema in the CMB have been extensively studied 
in the context of hotspots (peaks ) and coldspots (troughs) 
arising in Gaussian random fi elds l|Bond fc Efstathioulll987l: 
IVittorio fc Juszkiewicj [l987). Additional statistics related 
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to such local extrema, such as the Gaussian curvature 
and temperature-correlation function, have been investi- 
gated as a means to distinguish the geometry of the uni- 
verse and test the Gaussian hypothesis for the nature of 
the initial conditions jBarreiro et al.ll 19971 : iHeavens fc Shethl 



Il999l : IHeavens fc Gupta 2001 ). Such statistical techniques 
have subsequen tly been app l ied to seve ral datasets, and 
most notably bv lKogut et~aT1 (|l995l . ll99rj) with the COBE- 
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DMR data. Observations from the Wilkinson Microwave 
Anisotropy Probe ( WMAP) currently provide the most 
comprehensive, full-sky, high-resolution CMB measure- 
ments to-date, and studies on t he local extrema proper- 
ties thereof have been un dertaken (|Larson fc Wandeltlliooi . 
l2005l ; lToieiro et alj|2006h. 

More recently in lHou. Bandav fc Gorskl (|2009l ) . we ex- 
tensively analyzed the statistical properties of both the 
one- and two-point statistics of local extrema in the five- 
year WMAP temperature data. Such extrema are defined 
as those pixels whose temperature values are higher (max- 
ima) or lower (minima) than all of the adjacent pixels 
l|Wandelt. Hivon fc GorskH 1 199St ) . We considered only that 
part of the sky outside of the WMAP KQ75 mask and its 
subsequent partition in either Galactic or Ecliptic coordi- 
nates into northern and southern hemispheres (hereafter 
GN, GS, EN, ES). A frequentist comparison with the pre- 
dictions of a Gaussian isotropic cosmological model that 
adopted the best-fit parameters from the WMAP team was 
then made. The hypothesis test indicated a low-variance of 
both local maxima and minima in the Q-, V- and W-band 
data that was inconsistent with the Gaussian hypothesis at 
the 95% C.L. The two-point analysis showed that the ob- 
served temperature pair product at a given threshold v, 
^tt(9,v > 1,2), indicates a 3<r level 'suppression' on GN 
and EN, whereas £tt(# < 20°, v < oo) is suppressed on the 
full-sky and both northern hemisphere partitions. The latter 
is also the case for the point-point spatial pair-count function 
£pp(# < 20°, v > 1,2). Intriguingly, the statistics showed an 
^-dependence such that consistency with the Gaussian hy- 
pothesis was achieved once the first 5 or 10 best-fitting mul- 
tipoles were subtracted, implying that the anomalies may 
be connected to features of the large-scale multipoles. 

The local extrema anomalies therefore provide 
further evidence of a hemispherical asymmetry that 
was originally revealed using t he power sp e ctrum 
dHansen, Bandav fc Gorskil 12004 lEriksen et ail 12004 



lHansen et al . 
llEriksen et al.l 



2008), the N-point correlation functions 



2005) and Minkowski functionals (|Park C-Gl 



2004). This can be interpreted as a violati on of the cos- 
molog ical principle of isotropy. Theoretically, iGordon et al.l 
l|2005r ) proposed a mechanism of spontaneous isotropy 
breaking in which the long-wavelength modes of a medi- 
ating field couple non-linearly to the CMB perturbations. 
These fluctuations appear locally as a gradient and imprint 
a preferred direction on the sky. An implementation of a 
multiplicative modulation field of the intrinsic anisotropy 
with a single preferred direction was fitted t o the three-year 
WMAP observations in Eriks en et al. I i|2007l ) . The so-called 
'dipole modulation field' (dmf) was detected at a significant 
confidence level, and confirmed at even highe r signi ficance 
in the five-year WMAP data bv lHoftuft et all l|2009h . 

Another mech anism that violates rotational inva riance 
was proposed by lAckerman. Carroll fc Wise] <|2007l ). The 
ACW model picks out a preferred direction during cosmo- 
logical inflation to modify the power spectrum of primordial 
perturbations, P(k). Its imprint on the CMB can then be 
described by the covar iance matrix of spherical harm onic 
coefficients, (ap m a* ( i Tn i) . iGroeneboom fc Eriksenl (|2009l ) esti- 
mated the parameters of the ACW model in an extended 
CMB Gibbs sampling framework from the five-year WMAP 
data. The posterior distribution of the parameters indicated 



a convincing detection of isotropy violation at 3.8a signifi- 
cance in the W-band data. 

In this paper, we make a frequentist comparison of 
the dmf and ACW models with the five-year WMAP maps 
of temperature anisotropy using a large number of Monte- 
Carlo simulations. In particular, we consider the one- and 
two-point statistics of local extrema and evaluate whether 
the WMAP data are more consistent with these models 
as opposed to the standard Gaussian cosmological models 
against which various anomalies have been claimed. A rig- 
orous hypothesis test methodology is applied to establish the 
significance at which the models impact the previous conclu- 
sions regarding WMAP data and the violation of isotropy. In 
addition, the £tt and £pp correlation functions are utilised 
to further confirm the findings of the one-point analysis and 
to study the scale-dependence of the local extrema both in 
real and spherical harmonic space. 

It might be considered that the statistical significance of 
the local extrema anomalies alone is insufficient to warrant 
an analysis of the anisotropic models, in particular given 
the additional parameters required to specify them. How- 
ever, that they provide significantly better fits to the WMAP 
data than the isotropic model has been demonstrated by in- 
dependent Bayesian power spectrum analysis. It is clearly 
then quite legitimate to consider the local extrema statis- 
tics in terms of these improved models, and to determine 
whether the observed local extrema are also more consis- 
tent with such cosmological prescriptions. In fact, we will 
find that there is little evidence that the extrema anomalies 
are remedied by these models, and indeed in one case we 
show that the anomalous behaviour becomes more signifi- 
cant. 

This paper is organised as follows. In Section 12.11 we 
present an overview of the WMAP data and the instrumen- 
tal properties that must be involved in Gaussian simula- 
tions to enable an unbiased comparison with the real data. 
The two models and the algorithms for simulating them are 
briefly reviewed in Section 12.21 Section 12.31 prescribes the 
technique used for further data-processing and introduces 
the one- and two-point analysis methods. The results are 
reported in Section [3] before we present our conclusions in 
Section @] 



2 METHOD 

Both the dipole modulation and ACW models provide mech- 
anisms for the violation of cosmological isotropy in the con- 
text of Gaussian random fields. In a frequentist formalism, 
we construct simulations for each model, convolve with the 
appropriate instrumental beam transfer functions, then add 
noise according to the properties of the WMAP detectors. 
The statistical properties of the simulated data are then 
compared with the observed ones to verify the consistency 
of the models with the observed Universe. 



2.1 The WMAP data 

The WMAP instrument consists of 10 differencing as- 
semblies (DAs) cove ring frequencies from 23 to 94 GHz 
|Bennett et all 120031 ). The V and W-band foreground- 
reduced sky maps are used in our analysis, consistent with 
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the data sel ection for the WM AP five-year power spectrum 
estimation l|Nolta et al.l [2009) . The maps are available in 
the HEALPbH" 

pixelisation scheme with iV 8 i_e = 512 from 
the LAMBDA websitfl We combine data from DAs cov- 
ering the same frequency using uniform weighting over the 
sky and equal weight per DA. This allows a simple effective 
beam to be computed from the combination of beam trans- 
fer functions (be) for the DAs inv olved. In fact, the W MAP 
beams are not axisymmetric and Wehus et alj ([2009) have 
attempted to assess the effect of the combination of the true 
asymmetric beams with scanning strategy on the sky using 
reconstructed sky maps. They report deviations from the of- 
ficial bB at a level of wl.0% at I » 800. However, the addi- 
tional smoothing employed in this work render our statistics 
insensitive to such deviations, thus we continue to adopt the 
beam properties provided by the WMAP team. There are 
two noise components contributing to observations, white 
noise and residual 1/f noise following an effective destrip- 
ing in the map- making process. However, the 1/f noise can 
be considered negligible even on the l argest scales since th e 
5-year signal-to- noise level is so high l|Hinshaw et aJ.J[2007T ). 

The extended temperature analysis mask (KQ75 which 
leaves approximately 72% sky usable coverage) is applied 
to the data to eliminate regions contaminated by Galactic 
foreground and point sources. Hemispherical masks of the 
Galactic North and South (GN and GS), Ecliptic North and 
South (EN and ES) are also applied together with the KQ75 
mask for some parts of our analysis. In this paper, we will 
refer to results derived on the KQ75 sky-coverage as 'full- 
sky' results for convenience, to distinguish then from the 
corresponding results on masked hemispheres. 

2.2 The two models 

In this section, we briefly review the implementation of the 
two models discussed in our work. We use Td mf and Ta CW 
to identify the simulated realisations for each model, which 
are then combined with the beam profiles and WMAP-likc 
noise properties for each frequency band concerned in the 
analysis. The realisations are performed at a HEALPix res- 
olution parameter AT s ;de = 512 with maximum multipole, 
£ m ax = 1024. 

2.2.1 Dipole modulation 

We follow lEriksen et all (|2007f) and their implementation 
of a multiplicative dipole modulation field as proposed in 
iGordon et all i|2005l ). 

The CMB temperature field is defined as 

f(n) = T(n)[l + f(n)], (1) 

where T(n) = £\ aemYe m (n) is a statistically isotropic 
CMB temperature field with power spectrum Ce and f(n) 
is a dipole modulation field with amplitude much less than 

one but non-zero. 

In their low £ analysis of the WMAP data. lEriksen et ail 

1 http://hcalpix.jpl.nasa.gov/ 

2 http:/ /lambda. gsfc.nasa.gov/product/map/dr3/ 
maps_da_forered_r9_iqu_5yr_get.cfm 

3 http:/ /lambda. gsfc.nasa.gov/product/map/dr3/beam_info.cfm 



Table 1. The groups of parameters ado pted for the dmf si mula- 
tions. The (A, p) values are estimated bv lHoftuft et all ll2009h for 
different maximum values of multipole moment, ^ mo di involved 
during the estimation. The dipole axis p point towards (I, b) in 
Galactic coordinates. 

Band £ m od A p{l,b) 

0.119 - 0.034 
40 0.119 (224°, -22°) 

0.119 + 0.034 

0.080 - 0.021 
64 0.080 (232°, -22°) 

0.080 + 0.021 

0.070 - 0.019 
80 0.070 (235°, -17°) 

0.070 + 0.019 

0.074 - 0.021 
W 64 0.074 (232°, -22°) 

0.074 + 0.021 



(2007) have parameterised the dmf model using four free 
parameters, q (normalisation amplitude), n (tilt), A (am- 
plitude of dmf) and p (direction of dmf), together with a 
fiducial Gaussian power-spectrum, Cf d , which was chosen 
to be the WMAP best-fitting power spectrum. The power 
spectrum and modulation field are then written as 

Ge-,(q,n) = q ( jA Cf d and f(n;A,p) = An ■ p. (2) 

In fact, this formalism recognises that the best-fitting dmf 
may also impact other cosmological parameters, but that 
over the limited range in £ that the analysis is carried out, 
the practical effect can be represented by a normalisation 
q and spectral tilt n. The posterior distributions of the pa- 
rameters were then given by an optimal Bayesian analysis. 

In this work, we modulate Gaussian realisations of the 
CMB adopting the most recent paramet ers estimated from 
the WMAP5 data in lHoftuft et alj l|2009h . Their estimations 
seem to provide a robust value of p, whereas the magnitude 
A seems to be a function of (.mod, the maximum value of £ in- 
cluded in the analysis. We therefore simulate the dmf adopt- 
ing the best-fitting directions for the different £ nlo d ranges, 
and consider values for the amplitude A varying spanning 
the 68% confidence region for each range. These parameters 
sets are summarised in Table [T] Note that we apply these 
parameters over the full range of £ used in the simulations. 
Importantly, they find that q and n are consistent with unity 
and zero respectively (for the V-band, the actual values are 
q = 0.96 ± 0.05 and n = 0.06 ± 0.05), thus we can sim- 
ply assume the WMAP5 best-fitting 'lcdm+sz+lens' model 
power spectrum as our fiducial spectrum, and that the pres- 
ence of the dmf does not perturb the best-fit cosmological 
parameters. 

The simulated realisations T' drnl (n) are then given by 
a em = b t ^2f(n)Y em (n), (3) 

n 

TL(n) = ^a tm Y tm (n) + N(n). (4) 

These realisations are then convolved with the appropriate 
beam transfer function be, and noise realisations N(n) are 
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added. For each group of (A,p), we perform 2500 simula- 
tions for comparison with observations. 



2.2.2 The AC W model 

Inflationary theory gives a firm prediction that the primor- 
dial density pertur bations should be isotropic withou t any 
preferred direction. lAckerman. Carroll fc Wise] (1200 7T ) con- 
sidered the breaking of this rotational invariance by the ex- 
plicit introduction of a preferred direction, n' , during the 
inflationary era thus modifying the primordial power spec- 
trum P(k) to 



P'(k; g(k), n') = P(k) [1 + g(k)(k ■ n'/k) 2 



(•») 



where k is the vector with magnitude k. g(k) is a general 
function representing the rotationally non-invariant part, 
which ACW have argued can be approximated by a constant 
<7». The imprint of this preferred direction on the CMB is 
then given by 



Sl m ,l'm' — { a lm,at'm') — Ct^U'&n 



+ A 



£m,£' m' i 



(6) 



where Ce is the CMB angular power spectrum and we define 



9*& 



dkk 2 P(k)<di(k)<d e ,(k). (7) 



(,tm,£<m< are then geometric factors that describe coupling 
between modes, I to £'=£,£ + 2 and m to m' — m, m ± 
1, m±2. The integral term is a generalisation of the standard 
CMB power spectrum and may be computed using a mo di- 
fied version of CAMB (|Lewis. Challinor fc Lasenbvll2000h . 

The contribution of A can contribute significantly to 
the diagonal part of S which affects the total angular power 
spectrum and induces strong degeneracy between g* and the 
amplitude of the power spectrum of scalar perturbations, 
A s . However, it can be ensured that a given choice of g* 
mainly affects the anisotropic contribution by redefining the 
signal covariance matrix, as 



S, 



C£Sf e >S„ 



m.t' m' 



m.Sc ' m' 



l + g*/3 



(8) 



In this work, we assume that the cosmological param- 
eters are known and given by the WMA P5 best-fitting 
'lcdm +sz+lens' model power spectrum (IKomatsu et al.l 
2009). The two free parameters r emaining in the model, 

n '), have been estimated by iGroeneboom fc Eriksenl 
|2009l ) from the WMAP five -year data using Bayesian anal- 
ysis in an extended CMB Gibbs sampling framework. The 
estimation was carried out over different ^-ranges to investi- 
gate whether the preferred direction n' is dependent on the 
Grange selection. Indeed, a preferred direction and value 
for <j* for the ACW model was found for the V and W- 
bands at 2.5a and 3.8a confidence-level, respectively (Ta- 
ble 1 of their paper), although there is a 20° deviation in 
Galactic longitude between the two bands. For the range 
£ — [2, 400] and KQ85 sky-coverage, the g* value is reported 
as 0.10 ± 0.04 for the V-band and 0.15 ± 0.039 for W-band 
(95% confidence), while for the more conservative KQ75 
sky-coverage the V-band q* v alue is 0.10[— 0.100, 0.158]. In 
iHou. Bandav fc Gorskil (|2009l ). we show that the signal-to- 
noise-ratio (S/N) of the 5-year WMAP temperature data 
decreases to 1 at I = 399 (385) on V- (W-) band, and the 



Table 2. The groups of parameters used in our local extrema 
analysis for ACW model. The pa r amete rs follow the estimation 
results of lGroeneboom &: Eriksenl d2009t) for £-range [2,400] . 



V-band 
9* n'(l,b) 



W-band 
n'(l,b) 



-0.100 
-0.050 
0.050 
0.100 
0.158 



(130°, 10°) 



0.111 
0.130 
0.150 
0.170 
0.189 



(110°, 10°) 



0.158 (110°, 10°) 0.189 (130°, 10°) 



additional smoothing imposed during data-processing sup- 
presses information on this scale by a further factor of 0.1 
in spherical harmonic space. Therefore, we choose different 
g* values estimated from the ^-range [2,400] within the 95% 
confidence region, to carry out the frequentist comparison. 
The parameters used in this work ar e summarised in Ta ble [2] 
W e then follow the algorithm of lGroeneboom fc Eriksenl 
(2009) to get the simulated ACW realisations of the sky 
using 



Tl cw (n;g», n ) 



b.2~ 



E 



+ N(n), 



- tj, ( / 111 (g t , n')rn m Yi m (n) 



(9) 



where L^ m is the Cholesky decomposition of the sparse 
covariance matrix S, involving the best-fitting parameters 
(g*,n') and cosmological model, and r)i m is a set of Gaus- 
sian random numbers with zero mean and unit variance. The 

6 m0 -l 

factor 2 a is placed regarding the triangular Lt m and 
HEALPix data convention. Each simulation is convolved 
with the appropriate WMAP beams and noise realisations 
are added. For each group of parameters, we perform 5000 
simulations. 



2.3 Analysis and statistics 

In this section, we introduce the processing steps and statis- 
tics utilised for the local extrema comparison between the 
observed sky and simulations for our two selected models. 
Generally, the treatment follows closely t hat of sections 2.5 
and 2.6 of IHou. Bandav fc Gorskil (|2009h . but is briefly re- 
viewed here. 

The data and simulations are treated in an identical 
way. Masks are applied to the sky maps, then the best-fitting 
monopole and dipole components for the region outside the 
mask are subtracted - this is a standard procedure in CMB 
data analysis. Since our previous results on local extrema 
statistics suggest a close connection with large angular scale 
modes on the sky, we also consider cases where modes up 
to and including the quadrupole, £ ^ 5 or I ^ 10 are also 
subtracted before analysis. A modest smoothing is then ap- 
plied to the resulting sky maps where the angular scale of 
smoothing is determined by a signal-to-noise normalisation 
technique. In this work, we applied 43.485 arcmin and 45.064 
arcmin FWHM Gaussian beams to the masked sky maps of 
the V- and W-band respectively. In order to be conservative 
and avoid potential boundary effects in the analysis, the 
masks (for which all valid pixels are set to unity) are also 
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smoothed and only those pixels with values larger than 0.90 
are retained as valid. The local extrema are then determined 
using the hotspots program in the HEALPix package. 

First, the one-point statistics (number, mean, variance, 
skewness and kurtosis) of the local extrema temperature dis- 
tribution are computed for further hypo thesis test. The hy- 
pothe sis test methodology introduced bv lLarson fc WandeltJ 
(2005) is again adopted here. The probability of finding a 
simulated statistic that falls below the observed one is de- 
fined as p. As a two-sided test, we set the significance level 
a — 0.05 with the hypothesis 

H:p€ (a/2,1- a/2), (10) 

and perform the test twice to control the Type I error (/3 
- rejection of true hypothesis) and the Type II error (7 - 
acceptance of the false hypothesis) to be small. 

Secondly, the two-point correlation functions of temper- 
ature pairs, £tt, and spatial pair-counting, £pp, are com- 
puted to make a further x 2 analysis. During the two-point 
analysis, each realisation can be examined in a renormalised 
form, T(n) = ^(n)<T s k y , where <j s k y is the standard-deviation 
of the temperature field over the valid region of each reali- 
sation. We will continue to examine the correlation depen- 
dence on the thresholds u,e and the large-scale multipoles. 
For the £pp comput ation, we utilise the Hamilton estima- 
tor jHamijto^_^9^}_aiiduse_the same random sample as 
in iHou. Bandav fc Gorskl (|2009l ). Since the {£tt(6>)} and 
{£pp(#)} are asymmetrically distributed in a given angular 
range, the x 2 statistics are optimised by mapping the dis- 
tributions to Gaussian ones {s(8)}, 

Rank of observed map 1 f s _l* 2 , 

; ; ; — = / e 2 dt. (11) 

Total number of maps + 1 J_oo 

Then the probability of finding the number of simulations 
with a x 2 value below the observed one can be determined. 



3 RESULTS AND DISCUSSIONS 

In this paper, we will not attempt to describe all the one- and 
two-point results determined, but instead highlight some in- 
teresting cases that help to build the conclusions. 

Those results that are derived after subtraction of the 
low-£ multipoles in the range [0, £ T mv] are denoted by i Tmv 
for simplicity, where Irmv = 1, 2, 5 and 10. 

In the frequentist comparison, we determine the prob- 
ability p described in section [2~3l by counting the number of 
simulations with statistical values (either one-point statistics 
or the associated x 2 measure for the two-point analysis) be- 
low the observed one. For the one-point analysis, the proba- 
bilities are subjected to a hypothesis test twice, the first test 
sets (3 = a and any rejections are listed in the tables and 
marked by asterisks; the second test sets 7 = a, with the 
subsequent rejections that were accepted by the first test 
marked by question-marks. If we specify a = 0.05, those 
cases rejected by the first test indicate that the consistency 
of the observations with the model concerned is rejected at 
the 95% C.L., with an associated 5% probability of a Type 
I error (f3 = 0.05). For those cases accepted by the second 
test, it can then be asserted that such a consistency is ac- 
cepted at the 0.05 significance level, with an associated 5% 
probability of a Type II error (7 = 0.05). 



3.1 Dipole modulation results 

A subset of the results are shown in Table We only 
present the local extrema mean and variance statistics for 
Vp^'^SS) and W^= { 23 2 7 o 4 f_ 22 ^ here, since the number, 
skewness and kurtosis values are not revealing, and the re- 
sults from other groups of dmf parameters lead to similar 
conclusions. 

We first focus on the variance results from Table [3] For 
the most general case, l imv = 1, the probabilities for both 
the V and W-bands in the GN and EN are less-extreme than 
the unmodulated results as might be expected, although 
there are still rejections of the model on the full-sky and 
northern hemispheres. In addition, an improved agreement 
for l imv — 2 on EN and ES hemispheres indicates good con- 
sistency of the real data with the dmf model which again 
should not be surprising given its profile. The results on 
hemispheres implies that an increasing dmf amplitude may 
suppress the observed variance anomalies. However, after 
subtraction of the first 10 multipoles, new problems arise in 
the southern hemispheres, where the observed variance of 
local extrema is now much lower than the model prediction. 
This result is robust within the 68% confidence region of the 
dmf amplitudes tested here. It should also be noticed that, 
contrary to the cases on hemispheres, the full-sky variance 
is increasingly inconsistent with the model as the dmf am- 
plitude increases. This is suggestive that the dmf may not 
have the correct profile over the full-sky in order to reconcile 
the variance properties of the observed local extrema with 
the Gaussian anisotropic simulations. 

The observed mean statistics, which showed few anoma- 
lies when compared to the Gaussian isotropic model, are not 
consistent with the dmf, in particular for cases on the EN 
and ES. There are 3er-level inconsistencies in both the V and 
W-bands that demonstrate that the observed local extrema 
are not as extreme as the model expectations on the ES, but 
are too extreme on the EN, especially for £ lnlv — 5 and 10. 
Indeed, the discrepancies become more significant with an 
increasing dmf amplitude. It is noteworthy that the results 
for V^= 2 3 2 8 o f_ 22 o 2 ; and V^™±°$° (parameters fitted for 
(■mod = 64 and 80, respectively) show weaker evidence for a 
variance anomaly than found above, but the inconsistency 
of the mean values remains at 3cr C.L.. Therefore, we seem 
to have a 'paradoxical' situation where a weaker perturbed 
dmf is inadequate to explain the observed variance difficul- 
ties whereas the mean results and full-sky variances argue 
against a stronger dmf amplitude. 

The two-point results are more revealing. The full-sky 
results for both £tt(^<oo) and £pp(;/<oo) show few devi- 
ations from the non-modulated model and so cannot change 
the previous conclusions for these cases. For £tt(^>2), es- 
pecially with ^rmv = 5 and 10 as shown by the upper panel 
of Figure [1] the ^-frequencies of observation are increasing 
as the amplitude of the applied dmf increases. The relative 
spacing of the medians and the confidence interval bands 
shows that such a shifting is scale-dependent, decreasing on 
larger scales and providing evidence that the observation 
does not favor a full-sky dmf since, as found previously, the 
shape of the observed £tt(v > 2) fits the non-modulated 
median quite well. 

We then investigate £pp(f > 2) to find the impact on 
the spatial distribution of the local extrema by the dmf. The 
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Table 3. Frequencies of the extrema one-point statistics with lower values than the WMAP V and W-band five-year data. The dmf 
parameters are given as superscripts and subscripts of the frequency label. The local maxima (minima) results are denoted as 'max' 
('min'). NS, GN, GS, EN and ES correspond to full-sky, Galactic North, Galactic South, Ecliptic North and Ecliptic South sky-coverage 
outside of the KQ75 mask, respectively. The values rejected by the hypothesis test are marked by a * or ?. Values outside the 3<r 
confidence range are underlined. 
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lower panel of Figure [T] indicates that for realisations with 
the first 5 or 10 multipoles subtracted, the local extrema 
tend to be more clustered on scales less than 90° while more 
anti-clustered on larger scales, and this effect becomes more 
pronounced for a stronger dmf amplitude. However, the full- 
sky £pp(«/>2) medians in the equivalent unmodulated cases 
fit the observations quite well, and simply oscillate about the 
zero level on scales larger than 90° indicating an approxi- 
mately uniform spatial distribution of local extrema with 
v > 2. Again, the empirical properties of the dmf are not 
favoured by observations. 



The hemispherical results of £pp show similar behavior 
to the full-sky ones on scales less than 90° and so we focus on 
£tt results for the Ecliptic hemispheres here. As discussed 
in our previous work, on the GN and EN hemispheres with 
4 mv = 1 and 2, the observed ^tt(" < oo) is suppressed on 
scales less than 20° and £tt(v> 1,2) are disfavoured at the 
3<j-level in comparison with the Gaussian isotropic model. 
In contrast, the results for the GS and ES, as well as for 
correlation functions with £ IU1V — 5 and 10 over a variety 
of sky- coverages, show good consistency with the predicted 
medians. The results presented in Figure [2] here indicate an 
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Ecliptic North 



Ecliptic South 



180° 



Figure 1. The full-sky TT (upper panel) and PP (lower panel) 
correlation functions (180 bins) of local extrema for v > 2 
with £ rmv = 5 and 10. The red (blue) lined-dots correspond 
to the correlation function of local maxima (minima) observed 
in the WMAP V-band data. The light, middle and dark gray 
shaded bands show, respectively, the 68.26%, 95.44%, 99.74% con- 
fidence regions determined from 2500 MC simulations for dmf of 
V*jT(-224° -22°)' ant ^ bl ac k solid line shows the median, as 
well as the confidence regions and the median of +/ — la vari- 
ation of the dmf amplitude shown by orange/green lines (solid 
for the medians, dotted, dashed and dot-dashed for boundaries 
of 68.26%, 95.44% and 99.74% confidence regions, respectively). 
The probabilities of finding a simulation with corresponding \ 2 
values lower than the observed ones are marked for both maxima 
(red) and minima (blue). 
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Figure 2. The Ecliptic North and South £tt(^ < oo) (the up- 
per two panels for <? rmv = 2 and 10) and §tt(^ > 1) (the lower 
two panels for 4 mv = 1 and 10) with dmf V^H^q* ^. The 
nomenclature of the lined-dots, bands, lines and numbers follows 
the same style as Figure [JJ 



offset in the expected confidence regions between the north 
and south, and the discrepancy is enhanced with an increas- 
ing dmf amplitude. Considering the profile of the dmf, the 
model practically attempts to solve the suppression prob- 
lem on the northern hemispheres by imposing an effective 
north-south amplitude asymmetry that potentially improves 
the fit in that part of the sky. On scales less than 15°, the 
median of ^tt(^<oo) with £ ITav = 2 on EN does match the 
observations well. However, the model behaves poorly on the 
ES - the observations lie below the lower 3er-level for v < oo 
and below the 2<r-level for v>\. Moreover, after subtracting 
the first 5 or 10 multipoles, the north-south asymmetry is 
increased so far that the model cannot fit the observations, 
as quantified by the x 2 -frequencies and consistent with the 
behaviour of the one-point statistics. A simple interpreta- 
tion of these results would be that the profile of the dmf is 
incorrect, and that the amplitude of the effect may also be 
a function of £. 



3.2 ACW model results 

Generally speaking, the local extrema results in the ACW 
scenario show fe w differences from the corre sponding 
isotropic results in iHou. Bandav fc Gorskil (|2009l ) . For the 
northern hemispheres, the observed variances are still re- 
jected at 95% confidence level for i Tmv = 1 and 2, while other 
one-point statistics are quite consistent with the model ex- 
pectations, as judged by the hypothesis tests. Improvements 
in consistency occur again for <? rmv = 5 and 10. Nevertheless, 
we further investigate the variance results on the full-sky, 
GN and EN as a function of <?, (Figure [3} , and more so- 
phisticated acceptance/rejection boundaries are applied by 
setting different significance levels - a = 0.050, 0.025 and 
0.010. 

The p-values for the full-sky V-band variances exhibit 
an increasing trend g* value (note that we have added our 
previous results for isotropic (g» = 0.00) simulations to the 
plot), although such a trend is not apparent for the W-band. 
The frequencies of the £ IU1V = 1 minima (blue lines and dots 
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Figure 3. The frequencies (p) for which the local extrema vari- 
ances of 5000 ACW simulations fall below the WMAP5 V and 
W-band observations as a function of g t . Results are shown for 
the the full-sky (upper panel), GN (middle) and EN (lower). The 
red (blue) dots and lines correspond to the maxima (minima) 
frequencies for € rmv = 1, whilst the orange (green) values corre- 
spond to £ T mv = 2. The centres of the coloured circles correspond 
to the cases summarised in the last row of Table [2] Note that 
the red and blue circle in the lower-left panel overlap. The light, 
middle and dark grey-shaded bands fix the rejection regions in 
our two-sided hypothesis test with confidence level 95.0%, 97.5% 
and 99.0%, respectively, with the Type I error fj = a. The dotted, 
dashed and dot-dashed lines fix the lower-limit of the acceptance 
regions with a = 0.05, 0.025 and 0.01, respectively, with the Type 
II error 7 = a. A logarithmic coordinate is applied to the vertical 
axis to clearly show the acceptance/rejection boundaries. 



in Figure [3} are asserted to be rejected at 99% C.L., while 
not being 99% rejected (but still not accepted at the same 
level) when <;» increases to its 95% confidence upper limit, 
0.158. Despite the apparent increasing trend in the V-band, 
the £ rmv = 1 maxima (red lines and dots) variances are 
still not accepted at 99% C.L., while results of £ lmv = 2 are 
accepted within 95% C.L., also showing the increasing trend. 
This indicates that the ACW isotropy violation model shows 
mild improvement in the agreement with observations of the 
full-sky extrema variances. Nevertheless, the improvement 
for the cases with £ IU1V = 1 is inadequate to be accepted by 
the hypothesis test even at a significance level of 0.01. The 
results on GN and EN show similar features. The observed 
maxima variance on the W-band GN hemisphere increases 
with g* and is accepted at 0.01 significance level, as are the 
minima on the V-band EN hemisphere with £ rnlv = 2. There 
are no other acceptance/rejection changes as a function of 
9*- 

It is a natural concern as to the extent to which an incor- 
rectly aAa2t&A^v^06jmJ : h&jmeievred direction can affect re- 
sults. iGroeneboom fc Eriksenl (|2009T l found a 20° deviation 
of the best-fitting n' between V and W-band. Therefore, we 
assume that such a deviation represents the typical estima- 
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Figure 4. The full-sky and Galactic North £tt(^ > 1) (the upper 
two panels for £ rmv = 1 and 10) and £pp(v > 1) (the lower two 
panels for f rmv = 1 and 10) with 5000 ACW simulations on W- 
band, g* = 0.150 ± 0.039, n' = (130°, 10°). The nomenclature of 
the lined-dots, bands, lines and numbers follows the same style as 
Figure[T] Note that the medians and the boundaries of confidence 
region are more or less overlapped. 



tion error on the preferred direction and test its impact by, 
for example, fixing the maximum <?* estimated from the W- 
band but setting n' as that estimated from the V-band, as 
summarised in the last row of Tabled] However, as indicated 
by the coloured circles in Figure [3] no qualitative change is 
seen. 

Some of the two-point correlation results of ACW sim- 
ulations are highlighted in Figure [4] Even though the X*- 
frequencies decrease for larger values, the 3<r-level sup- 
pression of £tt(^> l,^rmv = 1,2) on the full-sky, GN and 
EN is still robust for the W-band, and the ACW model does 
not cause the ^-frequencies to improve sufficiently until the 
first 5 or 10 multipoles are subtracted. The confidence re- 
gions deviate little as varies within its 95% estimation 
error, constraining the observations almost at the same level 
as the isotropic model, for both £tt and £pp. It is thus con- 
vincingly indicated that the ACW model alone is unable to 
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reconcile the local extrema anomalies originally found in the 
context of an isotropic Gaussian model of CMB fluctuations. 



4 CONCLUSIONS 

In this paper, the statistical properties of local tempera- 
ture extrema in the 5-year WMAP data have been com- 
pared with two models that break rotational invariance - 
dipole modulation and the ACW model. Such a comparison 
was motivated by the fact that our previous analy sis of lo- 
cal extrema statistics l|Hou. Bandav fc~G 6rski 2009J) demon- 
strated an unlikely hemispherical asymmetry in the GN and 
EN as a consequence of an extremely low- variance compared 
to the expectations of a Gaussian isotropic scenario. More- 
over, the 5-year WMAP data indicate significant detections 
of these models on the basis of power spectrum analyses. 

We employ Gaussian MC simulations encoding the fea- 
tures of these two models to establish the statistical basis for 
testing whether the breaking of rotational invariance by the 
models can afford a satisfactory explanation of local extrema 
anomalies. As in previous work, both one and two-point 
statistics, as well as their dependence on large angular-scale 
modes have been studied. Both models are parameterised by 
a set of amplitudes and preferred dir ections as esta b lished 
by i ndependent Bayesian analyses in iHoftuft et al.l (|2009h 
and iGroeneboom fc Eriksenl (|2009T l . Our analysis has been 
carried out by sampling the former within their 95% confi- 
dence intervals, whilst imposing the fixed preferred direction 
found for each band, since the latter are estimated to be ro- 
bust for each model. The V and W-band data are considered 
here. 

In fact, the local extrema studies based on one and two- 
point statistics indicate that neither model provides a sat- 
isfactory solution to observed properties. In particular, the 
ACW model is probably the least interesting in this con- 
text - both the one-point analysis and two-point studies 
show similar features to the isotropic results, including their 
^-dependence. The dipolar modulation field, however, may 
itself be constrained by our analysis. 

Results determined from division of the data into hemi- 
spheres implies that a dmf with significant amplitude may 
suppress the observed variance anomalies. In particular, the 
4 mv = 2 values on EN and ES hemispheres indicates consis- 
tency of the data with the model. However, after subtraction 
of the first 10 multipoles, new problems arise in the southern 
hemispheres, where the observed variance of local extrema 
is now much lower than the model prediction. This is sug- 
gestive that the dmf may not have the correct profile over 
the full-sky in order to reconcile the variance properties of 
the observed local extrema with the Gaussian anisotropic 
simulations, and that the amplitude of the effect may also 
be a function of I. Moreover, the observed mean statistics, 
which showed few anomalies when compared to the Gaus- 
sian isotropic model, are not consistent and contradict the 
variance results by requiring a weaker amplitude. 

Our analysis indicates that neither a simple dipole 
modulated field model nor a rotational-invariance breaking 
model can satisfactorily explain the local extrema anomalies 
present in the WMAP data. On the one side, the rotational- 
invariance model affects the large-scale temperature ampli- 
tudes too little to significantly affect the local extrema statis- 



tics. On the other side, a modulation type model needs a 
more elaborate spatial structure than a simple dipole to fully 
fit the data. These issues remain interesting for further in- 
vestigation. 
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